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We describe two calculations involving P-wave mesons made of Wilson quarks: the strong coupling constant 
as in the presence of two flavors of light dynamical fermions and the mass and decay constant of the ai meson. 



1. INTRODUCTION 

P-wave meson spectroscopy has become in- 
creasingly interesting to the lattice community in 
the past few years, mainly for heavy quark sys- 
tems. The spectroscopy of heavy and light quark 
systems has its own intrinsic charm, and addition- 
ally the S-P mass splitting in quarkonium can be 
used to determine the QCD strong coupling con- 
stant through a technique pioneered by the Fer- 
milab group [1], where it is used to fix the lattice 
spacing and the strong coupling constant at a par- 
ticular q value is extracted from the expectation 
value of the plaquette. 

2. THE STRONG COUPLING CON- 
STANT 

We measure the strong coupling constant Ug 
from simulations using the HEMCGC /3 = 5.6 
configurations with two fiavors of light dynami- 
cal fermions [2], and from quenched simulations at 
l3 = 6.0 generated as part of the MILC collabora- 
tion program[3]. Our particles have valence Wil- 
son quarks, and so our study complements the 
work of the NRQCD collaboration[4], who used 
the same dynamical fermion simulations but non- 
relativistic heavy quarks. 

One begins by using the S-P mass difference of 
a heavy quark system to set the lattice spacing. 
In our case we have so far concentrated on the 
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= 442 MeV (1) 

where the ^Pi meson is the axial vector me- 
son created by the interpolating field i/'TsTsV'- 
We have also reconstructed masses of the full 
lattice multiplet corresponding to the P-wave 
system. This method allows us to compute 
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in addition to (1). Comparing lattice spacings 
obtained by fixing each of the above quantities 
gives a lattice spacing 1/a = 1900(50)(100) MeV, 
where the first error is statistical and the second 
is systematic. 

The coupling constant is defined through the 
plaquette in the so-called ay scheme. 
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Here \TrUp = 0.56500(2) so ay(3.41/a) = 
0.179(10), where the error is dominated entirely 
by 0{aY) effects. This is for two fiavors of dy- 
namical quarks. To convert to the physically 
interesting case (for charmonium) of three fia- 
vors, we repeat the calculation on the quenched 



I3 = 6.0 data set, where l/a = 2290(270) MeV 
and av(3.41/a) = 0.152(4), then run one data 
set to the same Q as the other and extrapolate 
linearly in l/ay from Nf = and 2 to Nf = 3: 

we find a^"^=^^(7.81 GeV) = 0.178(15). The 
largest contributions to the error come from two 
sources: the uncertainty in the lattice spacing and 
the Oiay) uncertainty in (4). 

Finally, to run to the Z mass (where numbers 
are conventionally compared) we convert from ay 
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and then run in Q using the formulas of Rodrigo 
and Santamaria[5]. Our result is 
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0.108(6). (6) 

A picture displaying our results with others at the 
Z mass is shown in Fig. 1. 

One puzzling feature of the calculation is the 
very different lattice spacings found by us using 
heavy Wilson quarks and the lattice spacing of 
the NRQCD group. The numbers are shown in 
Table 1, along with all other lattice spacings ex- 
tracted from this data set of which we are aware: 
in the table S and W label staggered and Wilson 
valence quarks, the "force" is from the string ten- 
sion [6], and the zero quark mass line is from an 
extrapolation when it is available. Clearly scal- 
ing violations are large in this data set. They are 
expected to be of order a for the hadron spec- 
troscopy, which presumably for this reason does 
not agree with experiment, and of order a^ for 
the "force" and NRQCD. 

3. THE ai MESON 

The ai meson is the isovector axial vector me- 
son made of nonstrange quarks. Long ago even 
the mass of the ai was controversial, but nowa- 
days the ai is most readily studied via tau decay. 
The coupling of an ai (of polarization fj) to the 
W-boson is parameterized[7] by the ai decay con- 
stant, fa^, by 

{0\i'Jij5lp\ai) = fa.Ci. (7) 

To calculate the mass and decay constant of the 
ai is a straightforward lattice exercise (provided. 
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Figure 1. Several measurements of the strong 
coupling constant at the Z-mass. The square 
marks our result. Diamonds represent other lat- 
tice calculations, and crosses are experimental re- 
sults. 
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Figure 2. The squares indicate measurements of 
j^iatt £qj. ^jjj.gg lig}it K values. The asterisk marks 
the extrapolation to Kcrit = 0.1610. 



Table 1 

Inverse lattice spacings (in MeV) from various observables from the dynamical fermion data sets. Statis- 
tical errors are 50 - 100 MeV. 
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0.025 1935 
0.010 2055 
0.0 2135 


2000 
2140 
2230 


1685 
1800 

1875 


1800 


1660 


2400 



of course that one can see a signal!) We merely 
carried the propagator calculations of Sec. 2 down 
to light quark mass. We used a combination of 
Coulomb-gauge Gaussian shell model sources and 
sinks and the axial current itself as a sink, and ex- 
tracted the mass and decay constants from corre- 
lated fits to the two propagators. We performed 
a jacknife extrapolation of the mass and decay 
constants into Kg, (see Fig. 2). 

We have so far only done this exercise for the 
l3 = 5.6 data sets. We have three values of light 
quark mass. We used our S-P mass splittings to 
set the lattice spacing. To convert the lattice de- 
cay constant to a continuum one we used the tad- 
pole improvement scheme of Lepage and Macken- 
zie[8]: 

(0|O-*(/. = l/a)\a,)j^={0\O'""{a)\a,)x 
/(m)a-2(am„j2(i + ^^^^ + 0(^2^) 

+Oia) (8) 

where f(m) = 1/4 (at k = Kg) converts the lattice 
field renormalization to the continuum, Aq = 0.31 
in tadpole improvement, and the factor a~^ con- 
verts the dimensionless lattice number to its con- 
tinuum result. The coupling constant is defined 
through the plaquette and is run down to a scale 

We did not see any dependence of our mea- 
sured quantities on quark mass, and so we present 
only averages here. We find nia-^ = 1270(70) MeV 
(expt. 1230(40)) and fa, = 0.31(3) GeV^ (expt. 
0.25(2)). Comparing with experiment is actu- 
ally rather nontrivial because of the large width 
of the ai; final state interactions are important. 
The "experimental" number for the decay con- 
stant comes from a phenomenological analysis by 
Isgur, Morningstar, and Reader[9], which uses ex- 
perimental input which is probably now obsoles- 
cent. 



4. CONCLUSIONS 

We plan to perform both of these measure- 
ments on the HEMCGC /3 = 5.3 dynamical Wil- 
son fermion configurations, and to complete a 
l3 = 6.0 quenched data set. These measurements 
will allow us to search for lattice spacing sys- 
tematics as well as possible effects of dynamical 
fermions. 
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